High-p T suppression and elliptic flow from radiative energy loss with realistic bulk 

medium expansion 
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We investigate nuclear suppression and elliptic flow in A + A reactions using Gyulassy-Levai-Vitev 
(GLV) radiative energy loss with the covariant transport MPC for bulk medium evolution. At both 
RHIC and LHC energies, we find that inclusion of realistic transverse expansion for the medium 
strongly suppresses elliptic flow at high pr compared to calculations with longitudinal Bjorken 
expansion only. We argue that this is a generic feature of GLV energy loss. Transverse expansion 
also enhances the high-pr suppression, while fluctuations in energy loss with the rescattering location 
of the jet parton in the medium lead to weaker suppression and smaller elliptic flow. However, unlike 
the strong reduction of elliptic flow with transverse expansion, these latter effects get nearly washed 
out once calculations are adjusted to reproduce Raa in central collisions. 

PACS numbers: 12.38.Mh, 24.85. +p, 25.75-q 



I. INTRODUCTION 

Understanding parton energy loss in ultrarelativistic 
heavy-ion reactions has been the focus of considerable 
recent theoretical effort. A variety of phenomenological 
approaches (e.g., [H-Q) formulate the problem determin- 
istically, in terms of a local energy loss rate dE/dL = 
—f(E(L),T(L),L) that depends on local temperature, 
position, and parton energy along the Eikonal (straight- 
line) parton trajectory. In the weak-coupling regime, 
more rigorous treatment is possible based on perturbative 
QCD^-Q- This includes quantum interference effects 
and also fluctuations, namely, energy loss for a given jet 
becomes a stochastic variable that is a function of the 
scattering and emission history of the jet. 

A critical step in computing heavy-ion observables 
from any energy loss model is spatial and temporal 
averaging over the bulk medium formed in the colli- 
sion. We employ here the Gyulassy-Levai-Vitev (GLV) 
framework Q in which a high-energy parton loses energy 
through gluon radiation induced by interactions with 
static Yukawa scatterers in the medium. It was natural 
for us to combine this approach with parton transport for 
the bulk evolution, namely Molnar's Parton CascadeM] 
(MPC). Selected early findings were highlighted in Q. 
Here we present a comprehensive set of results from a 
numerically improved calculation. 

Our approach is similar to recent work by Buzzatti 
and Gyulassy[l] or Horowitz [Toj], but with a few key dif- 
ferences. Unlike js}, we only focus here on light partons, 
and do not include multiple gluon radiation, elastic en- 
ergy loss, or energy loss fluctuations due to variations in 
radiated gluon momentum. However, as in Q, we do in- 
clude medium evolution with realistic 3D expansion, both 
longitudinal and transverse, which turns out to crucially 
influence elliptic flow and also affect the nuclear suppres- 
sion factor. We also study how differences between de- 
terministic (average) energy loss and stochastic energy 
loss, which fluctuates depending on where the jet parton 
interacts with the medium, impact observables. 



Recent assessment by the PHENIX Collaborationjl 1| 
highlighted the difficulty perturbative QCD parton en- 
ergy loss frameworks have with reproducing nuclear sup- 
pression at RHIC as a function of the angle with the 
reaction plane, Raa{4>)- GLV calculations were notice- 
ably absent from that analysis but our work here pro- 
vides a very similar cross-check of the GLV framework. 
The two observables we study, azimuthally averaged (tra- 
ditional) Raa arid elliptic flow i>2, carry essentially the 
same information (32j as Raa(<I>)- 



II. INGREDIENTS OF THE CALCULATION 

A. Energy loss framework 

We consider here the Gyulassy-Levai-Vitev (GLV) for- 
mulation of energy loss@, which provides the spectrum 
of gluons radiated by a jet parton as an expansion in 
the total number of scatterings n the jet and the radi- 
ated gluon experience with the medium as they travel 
through it. The dominant contribution is given by the 
leading n = 1 (single scattering) term 
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Here, x and k are the light-cone momentum fraction 
and transverse momentum of the radiated gluon, the jet 
energy is E, the original hard scattering is at z = 0, 
the jet is moving in the z direction and rescatters once 
at position z, /i(z) is the local Debye screening mass, 
X = J dz p a is the opacity for gluon jets (irrespectively 
of the jet Casimir Cr), and a = Qua 2 / (2p 2 ) is the 
(screened) total gg — > gg scattering cross section pro- 
vided the medium is made of gluons. The result was ob- 
tained for a medium of static Yukawa scatterers that are 
well-separated with mean distance d S> the radiated 
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gluon is assumed to be soft (x <C 1), and the calculation 
was done at tree level. Higher-order terms in number 
of scatterings have been systematically investigated (l^. 
up to at least n = 9, and were found to give modest 
corrections to the n = 1 result. 

We integrate the spectrum ([T]) numerically to obtain a 
momentum-averaged energy loss 
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for fixed z, i.e., retain energy loss fluctuations due to 
variations in z only. The probability for the scattering to 
occur at z is 
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Integrating over z as well yields the deterministic average 
energy loss 



(A£«) = J dzp(z)AE^( 
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For asymptotic jet energies, the average energy loss flU) 
is given by the GLV "pocket formula" 13 1 
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Though this remarkably compact expression (and its 
even simpler cousin without the ln[2i?/(/i 2 z)] factor) has 
motivated workers in the past, it is unreliable in prac- 
tice even for qualitative conclusions. Basic constraints 
that i) the radiated gluon energy does not exceed the 
energy of the parent jet, ii) momentum transfer in scat- 
tering off a thermal particle in the medium has on aver- 
age an upper bound given by the center-of-mass energy 
available, and iii) the radiated gluon energy must exceed 
the plasma frequency for the gluon to propagate in the 
medium, translate into kinematic bounds 



k <xE , q < V6ET , xE > fi 
and lead to (see App. |A"|) 



AE^(z) 



2Cna s 



E x l[b 



E 



-, e 



t(z) fi(z) 



(6) 



(7) 



where t(z) = 2E//j, 2 (z) is a characteristic formation 
time, and the function I is determined by the inte- 
gral (|A2jl . For small b — > and large E/fi — > oo, 
I(b,E/fi —> oo) sa — (7r/2)6 1nb, and one recovers ((3}. 
However, as shown in Fig. [I] the pocket formula does 
not account for smaller energy energy loss for less en- 
ergetic jets, and it also shows a peculiar maximum and 
turnover at higher z/t. 

Jet and medium parameters enter only through the 
arguments of I, which simplifies the calculation consid- 
erably because we can precalculate a 2D table and inter- 
polate later. Tabulating I is still time consuming because 



the integrand is oscillatory. To speed up the numerics, 
in H we used the approximate expression (|A6[) that in- 
volves one fewer integrals. Here we use the full result 
(|A2[) and confirm in Fig. [T] the accuracy of our earlier 
calculation (thin dotted lines). 
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FIG. 1: Function I(b, E / n) that governs GLV energy loss 
AE(z) in this work (see text). Simple analytic "pocket for- 
mula" (solid line) valid for asymptotic jet energies and small 
z/t is compared to the result with kinematic cutoffs (dashed 
lines) for various jet energies. Thin dotted lines are for an 
approximate calculation of 7 in our earlier work [j| . 

For very small b « 0, / can go slightly negative, im- 
plying a small energy gain instead of energy loss. In the 
calculation we cut these off and set I — whenever nega- 
tive contributions arise (see, e.g., the curve for Ej \x = 10 
below b < 0.1 in Fig.p. 

We emphasize that here we only consider radiative en- 
ergy loss as encoded in (fT]). There are significant addi- 
tional contributions from collisional energy loss[3|. Re- 
coil of scattering centers in scattering also enhances [l5| 
energy loss because it effectively changes the denomina- 
tor (q 2 + fi 2 ) 2 in (fT]) to q 2 (q 2 + fi 2 ). One also expects 
potentially large higher-order corrections in the strong 
coupling constant a s that are yet to be computed. In our 
approach all these effects get roughly incorporated into 
the parameter a s when we calibrate the calculation to a 
subset of observables through AE ~ a?. Equivalently, 
throughout this paper we apply (JTJ) with rescaled opaci- 
ties x Z\ with scaling factor Z — 3 (i.e., AE ~ Za 3 s ). 



B. Coupling to bulk medium 

Coupling energy loss to bulk medium evolution in- 
volves several open questions and assumptions (see, e.g., 
16] and references therein for a detailed discussion). 
First, the GLV result (fT]) was derived for a medium of 
static scattering centers whereas in a heavy-ion colli- 
sion the density changes rapidly with time. Unfortu- 
nately, the formulation has not been extended yet to 
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time-dependent background color fields. Therefore, as 
is customary in practice, for non-static media we rein- 
terpret p(z) in the GLV formula as the local density 
p(z, t = to + z) along the parton trajectory. 

Second, it is nontrivial to apply perturbative energy 
loss unless the medium is also perturbative. For exam- 
ple, questions arise whether the scatterers are quark- or 
gluon-like quasiparticles or some other field configura- 
tions, and whether the key parameter in ([1]), the local 
Debye mass /i ~ gT, is computed from the tempera- 
ture (~ T), local energy density (~ e 1 / 4 ), or some other 
thermodynamic quantity. For simplicity, we will con- 
sider here the medium to be a system of massless glu- 
ons that thermalize to a large degree by Bjorken proper 
time To = 0.6 fm. The local fi(z) is then determined 
along the straight-line jet path from the local tempera- 
ture, which we calculate using tabulated density evolu- 
tion n(xj_,?7 = 0,r) from the bulk dynamics model. For 
our gluon gas, n « 2T 3 and /i = gT s» 2T (similar ex- 
pressions have been used in [13]). We set aside the inter- 
esting question of quark chemical equilibration, namely, 
how quark/ antiquark abundances build up after the early 
gluon-dominated stage of the collision. At fixed parton 
density, quarks would significantly lower the opacity be- 
cause of their smaller Casimir (On the other hand, 
the Debye mass would only be affected by ~ f 0% based 
on perturbation theory [33l|.) 

Another open question is energy loss prior to the as- 
sumed thermalization time To. Reference [9] found GLV 
energy loss to be quite sensitive to whether density for 
t < To was taken to be constant, or rising linearly with t, 
or zero. However, differences in early density evolution 
were largely compensated by slightly different values of 
a s when the calculation was calibrated to experimental 
data. We assume here at early times t < tq a linear 
density build-up p(x.x,r),r) = Tp(x ± , rj, t )/t . 

One also has to address energy loss in the low- 
temperature (hadronic) regime, where the isolated 
Yukawa scattering center assumption in GLV clearly 
breaks down. In the calculation here we do not include 
hadrons, and simply set vanishing AE(z) = for jet path 
sections where fx < Aqcd, with cutoff Aqcd — 0.2 GeV. 

Though we expect that our qualitative conclusions are 
robust, it would be interesting to vary the assumptions 
above in the future. 



C. Jet initial conditions and fragmentation 

We take initial jet momentum distributions in 
p+p, Au+Au and Pb+Pb from leading-order (LO) 
perturbative QCD with one-loop running coupling 
a s(Q 2 ), using CTEQ5L parton distribution function 
parameterizations fl7j with Q 2 — Pr parton - Nuclear ef- 
fects such as shadowing are ignored but isospin (proton- 
neutron difference) is included. We consider jets pro- 
duced at midrapidity both in coordinate and momentum 
space, i.e., y = r/ = 0, where y and r\ are the coordi- 



nate and momentum rapidity 34]. The jet transverse mo- 
mentum distribution is generated uniformly in azimuth, 
while the transverse density distribution follows the bi- 
nary collision distribution for two Woods-Saxon distribu- 
tions pa{t)- Le., for collisions with impact parameter b 
we use the transverse profile 

1§7 = const x T -( x± + Ta ( x± ' (8) 

where Tyi(b) = J dzpA^V z 2 + b 2 ) is the nuclear thick- 
ness function. Woods-Saxon parameters for gold and 
lead nuclei were taken from HIJING 18]. To enhance 
statistics at high px, we generate uniform jet px in the 
intervals [2, 80] GeV for RHIC and [2, 300] GeV for LHC, 
and then appropriately reweight contributions to reflect 
the real jet spectrum. 

The density evolution of scatterers is calculated as dis- 
cussed in the next Section. After energy loss, jets are 
fragmented independently using LO BKK95 fragmen- 
tation function parameterizations [19j with scale factor 

Q 2 = PT,hadron> ancl we take n o = {^ + + k~)/2 for the 
neutral pion yield. This procedure reproduces high-p^ 
7To and charged particle spectra in p+p at RHIC and 
LHC with a modest if -factor Kmlo ~ 2.5 to account for 
higher-order contributions. 

D. Bulk medium evolution for r > to 

For medium evolution without transverse expansion, 
we use the initial conditions discussed later in this Section 
but keep the transverse profile "frozen", i.e., as in [9tll0j. 
the local density undergoes longitudinal Bjorken expan- 
sion p(x_l,t?,t) = p(xx,V> t o)to/t. 

For medium evolution with transverse expansion, we 
employ covariant transport theory in the same spirit 
as Refs. [2(J, Hl| to model the collective expansion of 
a causal, relativistic, low-viscosity bulk medium. We 
evolve a system of massless "gluons" with 2 — > 2 in- 
teractions via MPC 7] (Molnar's Parton Cascade), i.e., 
solve 

pffy/i = s(x, P i) + III [hU - hh) w 12 ^ 34 

234 

X 5 4 (pi+p 2 -p3~P4:) ■ (9) 

Here the integrals are shorthands for J. = J d 3 pi/ (2Ei), 
and fj = f(x,pj) is the phase space density. The source 
function S specifies the initial conditions, while W = 
(l/ir)s 2 do-/dt controls the local scattering rate. 

We use isotropic da/dt — 2a tot I s, for simplicity, with 
o~tot{ T ) = co(t/to) 2 / 3 growing with time so that the 
shear viscosity to entropy ratio rj/s stays approximately 
constant [U [22J. To generate substantial elliptic flow 
V2{pt ~ 3 GeV) ~ 0.25 in collisions with b = 8 fm im- 
pact parameter at RHIC and LHC, we set o-q = 8 mb 
and 4.5 mb respectively. 
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Initial conditions for Au+Au at ^/sWn = 200 GeV 
and Pb+Pb at ^snn = 2.76 TeV with impact param- 
eter 6 = 3 and 8 fm (about 0-10% and 30% centrality, 
respectively) are constructed as follows. As in [20[, we 
start at proper time tq = 0.6 fm with a locally ther- 
malized system at temperature T — 0.385 GeV. Because 
we are only interested in observables at midrapidity, we 
set up longitudinally boost invariant conditions in a wide 
coordinate rapidity window \r/\ < 5, with rapidity densi- 
ties dN{b)/drj oc N part (b) proportional to the number of 
participants ("wounded" nucleons) 

N part {b) = J d 2 x± [T A {x± + b) + T A (x± - b)] 

x (l - e - CT ««T A (x_ L ) N \ _( 10 ) 



Here the inelastic nucleon- nucleon cross section is ctmn — 
42 mb at ^/s^v = 200 GeV (RHIC) and 70 mb at 
y/SNN — 2.76 TeV (LHC), and we match the "gluon" 
multiplicity dN / dr\ to the observed charge particle dN/dy 
via setting dN(b = Q)/dr) = 1100 (Au+Au) and 2400 
(Pb+Pb) . For the transverse density profile, on the other 
hand, we use the same binary collision distributions as for 
jets. This choice is motivated by a classical Yang-Mills 
("color glass condensate") calculation 23] that found ec- 
centricities at early times to be much closer to binary 
collision eccentricities than to wounded nucleon eccen- 
tricities. 

The initial conditions above are more general than they 
may appear at first sight. Due to scalings[24j of the trans- 
port equation ©, a single transport solution contains the 
answer to a whole class of equivalent problems with scaled 
initial conditions, cross sections, and particle properties. 
For example, one can freely increase the initial density 
provided cross sections are reduced in inverse proportion 
to keep the mean free path the same. The initial temper- 
ature can also be rescaled together with particle masses 
(in fact for our massless quanta, To does not influence 
the density evolution at all). 

Though in principle the transport provides an ensem- 
ble of evolving scattering centers, here we solely use den- 
sity information to calculate energy loss. Technically this 
is very similar to employing ideal or viscous hydrody- 
namics for the medium evolution. We plan to explore 
different bulk dynamics models in the future. 



III. RESULTS 

Below we focus on two basic high-p^ observables for 
neutral pions at midrapidity, the nuclear suppression fac- 
tor Raa and the momentum anisotropy (elliptic flow) 
i'2 = (cos20) PT . Only energy loss is considered, i.e., con- 
tributions by the radiated gluons to the final spectrum 
and feedback on the bulk medium due to the jet were ig- 
nored. This is a good approximation at sufficiently high 
Pt- 



We consider four scenarios based on i) whether the 
medium is only undergoing Bjorken expansion ("ID" as 
in d,[I3|) or transverse expansion as well ("3D"); and ii) 
whether average energy loss (AE) is used or the stochas- 
tic AE{z). 

Before we turn to results for RHIC and LHC, it is illus- 
trative to compare these scenarios for gluon jets of fixed 
initial energy in Au + Au at RHIC with 6 = 8 fm. Every 
aspect of the calculation (such as production points with 
binary collision profile, uniform angles in azimuth, etc) is 
as discussed above except that the gluon jets all start out 
with energy Eq = 20 GeV. Figure shows the final jet 
energy distribution after energy loss. The first obvious 
feature is that in both ID and 3D scenarios, stochastic 
energy loss gives a broader energy loss distribution than 
using the average (E) . This is natural since averaging re- 
duces fluctuations in general. Given that parton spectra 
are convex at high px, more fluctuation with stochastic 
energy loss should result in weaker nuclear suppression 
(higher Raa), for the same a s . 
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FIG. 2: Final energy distribution for gluon jets of fixed ini- 
tial energy 20 GeV in four scenarios (see text) for Au + Au at 
^snn = 200 GeV with 6 = 8 fm, using averaged or stochastic 
energy (solid vs dashed) with or without transverse expansion 
for the medium (lines without symbols vs lines with symbols). 
All distributions are normalized to unit area. In the 3D case, 
the medium evolution was computed using the parton trans- 
port MPC0- 

Another noticeable feature in Fig. [5] is that the distri- 
bution is narrower with transverse expansion (3D) than 
with longitudinal expansion only (ID). This is a result 
of the interplay between energy loss and transverse ex- 
pansion of the density profile. In the "frozen" ID case, 
largest energy loss occurs for jets that start spatially out 
of plane and move through the center of the collision zone 
along the long axis of the almond-shaped transverse pro- 
file. With transverse expansion, however, the situation 
is the opposite. Largest energy loss occurs for jets that 
start in-plane and move along the initially shorter axis of 
the almond. Because the medium expands faster in the 
in-plane direction, by the time these jets cross the center, 
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the spatial eccentricity largely disappears (near the colli- 
sion center the elongation is actually in-plane already but 
the tail of the density distribution away from the center 
still keeps traditional eccentricity e = (y 2 ~ x 2 )/(x 2 + y 2 ) 
positive). The low-energy side of the final jet energy dis- 
tribution in Fig. [2] cuts off more quickly in the 3D case 
because the probability for hard scattering (forming a 
jet) at the edge of the collision zone drops much sharper 
in-plane than out-of-plane. 

Rapid transformation of the collision zone towards a 
more axially symmetric shape also reduces the energy 
loss difference between jets crossing through the collision 
center in-plane vs out-of plane (what in-plane jets miss 
in opacity while approaching the center they make up for 
while leaving it, and vice versa), which explains why the 
energy loss distribution is much flatter in the transversely 
expanding case. 



A. Nuclear suppression in Au+Au at RHIC 

Figure [3] shows our results for neutral pion Raa 
at RHIC, for fixed a s = 0.29 for all four scenarios. 
With stochastic energy loss the suppression is noticeably 
weaker, as expected for "upward curving" parton spec- 
tra at high pt- Quite interestingly, realistic transverse 
expansion significantly enhances jet quenching. This is 
a generic feature of GLV energy loss coming from the 
(1 — cos) interference term in (TT]). Scatterings at large z 
induce larger energy loss (cf. Fig. [IJ , and with a trans- 
versely expanding density profile there is higher chance 
to scatter further away from the production point than 
in the transversely frozen case. 

Unfortunately, without precise control over a s , Raa 
alone cannot differentiate between our four scenarios. As 
shown in Fig. [4] after a slight tuning of a s to reproduce 
the suppression in central collisions, differences in Raa 
largely disappear. We need ~ 10 — 20% higher a s with 
fluctuating energy loss than with the path-averaged one, 
while ~ 10 — 15% lower a s with transverse expansion 
than without it. 



B. Elliptic flow in Au+Au at RHIC 

Luckily, p^-differential elliptic flow turns out to be a 
much more sensitive probe of the medium evolution and 
energy loss treatment. Figure [5] shows our results for 
neutral pion V2(pt) in Au + Au at y / 5jv r /v = 200 GeV at 
RHIC, for fixed a s — 0.29 in all four scenarios. We find 
that, especially in more peripheral b — 8 fm collisions, 
stochastic energy loss gives smaller elliptic flow both with 
and without transverse expansion than average energy 
loss. This is in line with the weaker nuclear suppression 
for the stochastic case shown in Fig. [3] 

On the other hand, transverse expansion reduces ellip- 
tic flow, and at the same time also gives smaller Raa 
as we have seen in the previous Section. This may seem 
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FIG. 3: Neutral pion Raa in Au + Au at y/s^N = 200 GeV at 
RHIC with impact parameters b = 3 f m and 8 fm (lower pair 
of lines vs upper pair of lines, respectively, in both panels) 
calculated in four different scenarios with GLV energy loss 
(see text). Left panel: transversely frozen density profiles. 
Right panel: realistic transverse expansion modeled with par- 
ton transport MPC@|. Dashed curves were obtained with 
fluctuating energy loss, while solid curves with averaged en- 
ergy loss along the jet path. Data[25l] from PHENIX (boxes) 
are also shown for comparable centralities 0-10% and w 30% 
to guide the eye. 
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FIG. 4: Same as Fig. [3]but with a 3 tuned in all four scenarios 
to reproduce Raa in nearly central (6 = 3 fm) collisions. 



counter- intuitive at first, but it is also a feature of GLV 
energy loss. Because the interference term in (JTJ biases 
against early scattering in the medium, a jet parton loses 
more energy if it scatters further away from the produc- 
tion point, which also means the scattering is later in 
time. Compared to the transversely "frozen" ID sce- 
nario, transverse expansion gives at later times higher 
densities away from the center, so there is more quench- 
ing, but at the same time the medium also rapidly be- 
comes more azimuthally symmetric, so the elliptic flow 
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response is weaker. 
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FIG. 5: Neutral pion V2 in Au + Au at ^/snn = 200 GeV at 
RHIC with impact parameters b = 3 f m and 8 fm (lower pair 
of lines vs upper pair of lines, respectively, in both panels) 
calculated in four different scenarios with GLV energy loss 
(see text). Left panel: transversely frozen density profiles. 
Right panel: realistic transverse expansion modeled with par- 
ton transport MPC0]. Dashed curves were obtained with 
fluctuating energy loss, while solid curves with averaged en- 
ergy loss along the jet path. Dataj2(J from PHENIX (boxes) 
are also shown for comparable centralities 0-10% and w 30% 
to guide the eye. 

The reduction of V2 in the 3D case is manifest even af- 
ter a s is tuned to reproduce Raa in central collision, as 
shown in Fig. |6l For both averaged energy and stochas- 
tic energy loss, v-i is reduced by 40% (nearly half!) 
at high pt ~ 10 — 15 GeV due to transverse expansion. 
Similar reduction is present at LHC energies as we dis- 
cuss in Section UlI Dl In the more realistic "3D" scenario 
with both transverse and longitudinal expansion, both 
averaged and stochastic energy loss give about the same 
V2(pt) a t high p? for both centralities we studied. 



C. Nuclear suppression in Pb+Pb at the LHC 
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FIG. 6: Same as Fig. [5]but with a 3 tuned in all four scenarios 
to reproduce Raa in nearly central (6 = 3 fm) collisions. 



cases studied stochastic energy loss gives somewhat flat- 
ter Raa than a computation using the average energy 
loss. From our fixed-coupling calculation, Raa curves 
rise with px slower than the trend in the data. Incor- 
poration of running a s (Q 2 ), as was done, e.g., in [Hj], 
should yield better agreement. 
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Now we turn to results for Pb+Pb at LHC energies. 
For the nuclear suppression, we find the same generic fea- 
tures as in Au+Au at RHIC: i) medium evolution with 
transverse expansion gives stronger suppression (smaller 
Raa) than the transversely "frozen" scenario, and ii) 
stochastic energy loss gives less suppression than the 
path-averaged (AE 1 ). 

Figure [7] shows our results for neutral pion Raa in 
Pb+Pb at y/sWw = 2.76 TeV at the LHC, after a s has 
been adjusted to reproduce the suppression measured 
around pr ~ 8 GeV. At LHC energies we only computed 
the transport evolution for b = 8 fm collisions, therefore 
in the 3D scenarios we adjust a s to Raa measurements 
for « 30% centrality(6 = 8 fm) instead of — 10% cen- 
tral collisions (6 « 3 fm). We generally find that in all 



FIG. 7: Neutral pion Raa in Pb + Pb at ^/sWn = 2.76 TeV at 
the LHC with impact parameters b = 3 f m (left panel only) 
and 8 fm (both panels) calculated in four different scenar- 
ios with GLV energy loss (see text). Left panel: transversely 
frozen density profiles. Right panel: realistic transverse ex- 
pansion modeled with parton transport MPC0]. Dashed 
curves were obtained with fluctuating energy loss, while solid 
curves with averaged energy loss along the jet path. Charged 
hadron data[27| from CMS (boxes) are also shown, scaled to 
comparable centralities 0-10% and « 30%, to guide the eye. 

We need somewhat larger a s values to match the ob- 
served suppression at RHIC than to reproduce suppres- 
sion at the LHC (Fig.@]vs Fig. [7]), i.e., the medium at the 
LHC appears somewhat less opaque to jets than extrap- 
olation based on RHIC data would suggest. Our results 
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thus qualitatively reinforce the findings in [29], although 
we see about twice as large an effect. In our transversely 
static (ID) scenarios the ratio of LHC to RHIC effective 
couplings is about a^ HC / af HIG ps 0.8, while with trans- 
versely expanding medium the ratio is somewhat lower, 
a^ HC /af HIC ~ 0.75. 



D. Elliptic flow in Pb+Pb at the LHC 



Finally we demonstrate that with GLV energy loss the 
striking reduction of high-p^ elliptic flow seen at RHIC 
energies for transversely expanding media also occurs 
at LHC energies. Figure [5] shows our results for neu- 
tral pion V2{pt) in Pb + Pb at ^/snn — 2.76 TeV at 
the LHC, again after a s has been adjusted to reproduce 
Raa{pt ~ 8 GeV). We only find a very small reduction 
of elliptic flow due to stochastic energy loss compared 
to using the average (AE) along the jet path. The dif- 
ference between transversely "frozen" evolution and full 
3D expansion, however, is even larger than in Fig. [5] for 
RHIC collisions. By px — 5 GeV, i>2 is reduced by more 
than half, and above px = 15 GeV by more than two- 
thirds, irrespectively of whether we use average energy 
loss or the stochastic one. 



3D, stochastic 

3D, averaged 
ALICE h 



o 




15 25 35 
P T [GeV] 



15 25 35 
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FIG. 8: Neutral pion vi (p T ) in Pb + Pb at y/stm = 2.76 TeV 
at the LHC with impact parameters b — 3 fin (left panel only) 
and 8 fm (both panels) calculated in four different scenarios 
with GLV energy loss (see text). Left panel: transversely 
frozen density profiles. Right panel: realistic transverse ex- 
pansion modeled with parton transport MPC[3]. Dashed 
curves were obtained with fluctuating energy loss, while solid 
curves with ave rag ed energy loss along the jet path. Charged 
hadron data[30l.l3"H from CMS (boxes) and ALICE (triangles) 
are also shown, scaled to comparable centralities 0-10% and 
« 30%, to guide the eye. 



IV. CONCLUSIONS 

In this work we investigated nuclear suppression and el- 
liptic flow in A + A reactions at RHIC and the LHC using 
Gyulassy-Levai-Vitev (GLV) energy loss with the covari- 
ant transport model MPC[7j for the bulk medium. We 
compared calculations with transversely "frozen" density 
profiles (as in 0, [lj| ) to calculations with realistic trans- 
verse expansion, and also studied the difference between 
path-averaged (deterministic) energy loss and stochastic 
energy loss that fluctuates depending on the location of 
the scattering center that interacts with the jet parton. 

Our most striking finding is that, at both RHIC and 
LHC, realistic transverse expansion strongly suppresses 
elliptic flow at high px compared to calculations with 
longitudinal Bjorken expansion only. We argue that this 
is a generic feature of GLV energy loss, coming from in- 
terference bias against early rescattering in the medium. 

We also find that transverse expansion enhances high- 
Pt suppression, while energy loss fluctuations in the 
stochastic case lead to weaker suppression and smaller 
elliptic flow. However, unlike the strong reduction of el- 
liptic flow with transverse expansion, these latter effects 
nearly disappear once calculations are adjusted to repro- 
duce Raa in central collisions. 

Though our calculation lacks a few relevant effects such 
as elastic energy loss and multi-gluon emission, the re- 
sults do suggest that GLV energy loss will have difficulty 
with the simultaneous description of nuclear suppression 
and elliptic flow at RHIC and the LHC. In view of the 
Raa{&) challenge posed for perturbative QCD energy 
loss in a recent compilation [Tl| by PHENIX, we believe 
it is imperative to test our findings with more full-fledged 
GLV implementations foL Hoj. and also other bulk evolu- 
tion models (such as hydrodynamics). 
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Appendix A: The energy loss integral / 

Here we discuss how Eq. |(7J) is obtained. Substituting 
([1]) and ([6]) into ([2]) leaves one with five integrals to do. 
Axial symmetry cPqcPkdx = (tt/2) d(q 2 ) d(k 2 ) dipdx re- 
duces these to four, and with dimensionless variables 



k 



E 
/'- 



(Al) 



8 



and u = (k — q) 2 / fj 2 = ft 2 + £ 2 — 2k£ cos y?, we have 

3e 7T e 1 



xe 7T 



km - /«« 2 >/^/«« / 



x mi „ 

x I 1 — COS 

a; 



2£ cost^ 



where the lower limit on the x integral is x m i n = 
max(l,ft)/e. The innermost integral is doable with the 
help of 



, ol a _,.q 

ax cos — = x cos — h ev Si— 

x xx 



(A3) 



where Si is the sine integral function. The remaining 
three integrals can be evaluated numerically (at some 
expense because the integrand is oscillatory). 

For a numerically easier approximate result, one may 
momentarily ignore the upper bound © on q. With new 
variable q — > Q = (k — q)//i, Q is then unrestricted 
and we have, again in polar coordinates d 2 Qd 2 krfa; = 
{it/2)d{Q 2 )d(k 2 )d<j)dx, 



x / dn 



2(k — Qcos( 



7T (k 2 + u — 2kQ COS(f) + 1) 







(A4) 



(A2) The <p and k integrals evaluate to 



g{u,xe) = 



1 + u ^/ x 2 e 2 + 2xe(l - u) + (1 + uf 



(A5) 



Restoring now an approximate ^-averaged upper limit 
on Q 2 , Q 2 = u ~ k 2 + £ 2 < x 2 e 2 + 3e < Axe 2 gives the 
simpler 2D integral 



1 4xe 

h D (p,e) = J dx J 



ub . 

cos — gtu, xe) , (A6) 
x 



l/e 



I(b, e) -> JdxJ—^l-cos — 

1/6 



which approximates the full result (| A2|) well in practice 
provided E/fi> 10 (cf. Fig. [TJ. 
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